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^ ' Abstract 

> : 

^ ■ This report consists of three parts. In Chapter 1 we present a brief description of the 

LEPTOP approach for the calculation of the radiative corrections in the Minimal Stan- 
dard Model. The approach is based on the one-loop approximation with respect to the 
genuinely electroweak interaction for which simple explicit analytical formulae are valid. 
Starting from G^,mz and a = a(m 2 z ) as the main input parameters, we consider all 
observables of the Z boson decays and the mass of the W boson on the same footing. 
The dressing with gluonic corrections is performed by using the results of calculations 
published by other authors. 

In Chapter 2 theoretical uncertainties inherent to our approach and hence to the LEPTOP 
computer program are analyzed. 

In Chapter 3 we describe the Fortran code of LEPTOP. This code can be obtained on 
request from rozanov@afsmail.cern.ch. 

This document can be accessed on http : // cppm . in2p3 . f r/leptop/intro_leptop . html 

by WWW users. 



Chapter 1: LEPTOP BASIC EQUATIONS 
1.1 Introduction 

The program LEPTOP calculates radiative corrections to the observables of Z boson 
decay and of the W boson mass, mw, in the Standard Model |L|. It is based on a set of 
explicit formulae for these observables, written in the papers - 0. 

Our aim is the study of the genuinely electroweak (e.-w.) interaction by comparing the- 
oretical results with experimental data. In LEPTOP we separate the e.-w. interaction 
from the well known pure electromagnetic effects. Such separation is well defined in 
the one loop e.-w. approximation. Then we include all gluonic corrections (internal for 
hadron-free observables and both internal and external for the Z boson hadronic decays) 
known in the literature. The main result of our study (refs. is that within la 

- 2a accuracy experimental data are described by the e.-w. Born approximation which 
is universal for any modification of the Minimal Standard Model. The one-loop e.-w. 
corrections (that really specify the e.-w. theory) are barely visible at the level of LEP1, 
SLD, UA2 and CDF accuracy. We do not need to consider two-loop e.-w. corrections, 



except for enhanced gaugeless top-quark contribution [|T0||- [[13 

In LEPTOP we confine ourselves by considering effects larger than 10 -5 . 

Renormalization in LEPTOP . Loop calculations depend on parameter e = 4 — D in di- 
mensional regularization (D is dimension of space-time), on the 't Hooft's scale parameter 
fi, on two bare gauge coupling constants, vacuum expectation value of the Higgs field and 
on the bare masses of quarks, leptons and the Higgs boson.0 The dependence on e and 
H disappears if one rewrites the result of loop calculations in terms of renormalized cou- 
pling constants, physical masses, etc., or in terms of corresponding number of some other 
independent physical observables. 

Input parameters. In LEPTOP we present the results of the loop calculations in terms of 



three most accurately known electroweak parameters [13], | ]15| |: 

m z = 91.1887(44) GeV (1.1) 

Gf, = 1.16639(2) • 1CT 5 GeV~ 2 (1.2) 

a = a(m%) = 1/128.87(12) (1.3) 

and three free unknown parameters: physical top quark mass rrit (on mass shell), higgs 
mass m H and strong coupling constant a s at scale m z : a s = a s (q 2 = m 2 z )jjg] the latter 
enters through gluonic corrections. There are tiny effects of nonzero 6-quark mass rrib, 
r-lepton mass m T and c-quark mass m c . All other quarks and leptons can be considered 
as massless. 

Note that G M , a and a s are not directly observable quantities: 

1) Gfj, is a four-fermion coupling constant of the muon decay corrected by purely electro- 
magnetic loops. 



1 In general they also depend on Kobayashi-Maskawa angles but in the limit mj 1 /m 2 z 1 for all quarks 
q =/= t this dependence vanishes for Z boson decay and mw- 
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2) Electromagnetic coupling constant a at scale q 2 = m| includes running from the value 
a to a due to the three lepton and five "light" quark loops calculated by formulae: 

m| f°° ds . . _ 

da = — r — a e+e -(s) (1.5) 

4w 2 a Jthr m z — s 

where a is fine structure constant and u e + e - is the cross section of e + e~-annihilation, 
through one virtual photon, into leptons and light hadrons: 

5a = 5a t + 5a h = 0.0314 + 0.0282(9) . (1.6) 
The leptonic contribution 

5a, = ^S ; (ln - f ) = 0.0314(0) (1.7) 
Sir mf 3 

is known to high accuracy. The hadronic contribution, Sa^, is calculated by using eq. 
( |1.5|) with the experimental cross-section a e + e -^h(s) in the dispersion integral below s = 
(40GeV) 2 and parton result above s = (40GeV) 2 (see ref. fT5j ). 



3) The strong coupling constant a s at scale q 2 = m 2 z is defined as renormalized coupling 
in the MS scheme of subtraction. 

By using a instead of a we automatically take into account the fact that the running 
weak coupling constants a>w{.(f) and az(q 2 ) do not run actually in the region \q 2 \ < m 2 z , 
while a(q 2 ) does run. So a is relevant to electroweak corrections (not a in spite of its 
extremely high accuracy) @, 0. 

The outline of Chapter 1 is as follows: In sect. 2 we introduce the general notations and 
phenomenological expressions for observables including the final state interaction due to 
photon and gluon exchange up to the terms of the order of 10~ 5 . In sect. 3 we consider 
the a-Born approximation. One e.-w. loop corrections are presented in sect. 4 for hadron- 
free observables (sect. 4.1) and for Z boson decay into hadrons (sect. 4.2). Chapter 2 is 
devoted to the accuracy of the LEPTOP. In sect. 5 we make general remarks concerning 
various sources of uncertainties. In sect. 6 we consider uncertainties in the vector boson 
self-energies. Sect. 7 is devoted to the uncertainties in the hadronic Z decays. The 
procedure of estimating the total theoretical accuracy of LEPTOP is described in sect. 
8. Appendix A presents the flowchart of LEPTOP. Appendix B contains the full list of 
our papers in which the LEPTOP approach was developed. 

1.2 Phenomenology and notations 

It is useful to define the electroweak angle 9 (sin 8 = s, cos 9 = c) in terms of three basic 
parameters , m z and a: 

2„2 1 o 
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(for earlier references see [16]). Solving eq. ( PL.8Q with experimental values eqs. (1.1)— (1.3) 
we get 

s 2 = 0.23117(33) 

c = 0.87683(19) (1.9) 

In LEPTOP the amplitude for Z boson decay into fermion-antifermion pair // is written 
in one e.-w. loop approximation in the following form 

f 2 = W^m! = 0.54866(8) (1.10) 

In this amplitude we neglect such structures as the weak magnetic moment induced by 
e.-w. interaction. These structures are of the order of (a/7rs 2 )(mf/mz) 2 , where mj is 
the fermion mass, that is far beyond the accuracy of the present and future experiments; 
even for 6-quark it is ~ 10~ 5 . 

By definition the constants gy.A do not include the contribution from the interaction in 
the final state due to the exchange of gluons (for quarks) and photons (for quarks and 
leptons). The final state interaction includes also the processes of photon and gluon 
emission. Corresponding corrections have nothing to do with electroweak corrections and 
can be written separately. 
Leptonic widths 

For the decay into charged leptons 11 we explicitly take into account the final state QED 
interaction factor 



4To[(^) a (i + ^a) + {g l A f{i + - 6^)] (l.n) 



where 



r = ^fa G i* m z = 82.945(12)MeV . (1.12) 



Second order QED corrections (-) 2 ~ 10 6 are neglected. The term proportional to 
leptonic mass squared in eq. Ql.llj) is also negligible for I = e, \i (so we put m e = m M = 0) 



and is barely visible only for I = t (m 2 /m 2 z = 3.8-10 4 ). 
For the neutrino decay 

= 8<^r , 

g v = gvu = gAv ■ (1.13) 

Hadronic widths 

For the decay into light quarks q = u,d, s we neglect small effects due to nonzero quark 
mass (i.e. we put m u = = m s = 0) and take into account the final state gluon 
exchange up to the third order |T7|], JT8], |T9|, the final state one photon exchange 



and the interference of the photon and gluon exchange |21J . These corrections are slightly 
different for vector and axial channels. 
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For the decay into quarks we have 

F q = T(Z - qq) = 12[g 2 Aq R Aq + g 2 Vq R Vq ]V Q (1.14) 

where factors Ra,v are due to the final state interaction. (In our previous papers we used 
letter G instead of R for these factors). According to |18|], |nj, |20 



71 4 y 7T 4 y 7T 7T 



+ [1.409 + (0.065 + 0.015 lnt)-](—) 2 -12.77(—) 3 + 12— f— 5„ m (1.15) 

t 7T 7T 7T 

= fly, - (2T 3g )[/ 2 (t)(^) 2 + J 3 (t)(-) 3 ]- 

7T 7T 

-12— f— 6 vm - 6^6 am - 10— (— ) <5 am 1.16 



where m q is running quark mass (see below) 

■ Oils ._ . -,OL 



5 vm = 1 + 8.7(— ) + 45.15(— Y, (1.17) 

7T 7T 



*L = 1 + 3-67(-) + (11-29 - lnt)(^) 2 , (1.18) 

7T 7T 



Q 8 \nt , 

/,(«) = -3.083 -ln i + ^ + °™ , (1.20) 

t V 



J 3 (t) = -15.988 - 3.722 hit + 1.9171n^ , (1.21) 

t = m 2 t /m 2 z . 

The terms of the order of (— ) 3 due to diagrams with three gluons in the intermediate 
state have been calculated for Ry q in ref.|2(J. They have small numerical coefficients so 



that the net effect is of the order of 10 5 and we have omitted these terms in expansion 



For Z — > bb decay the nonzero 6-quark mass is not negligible and produces contribution 
of the order of 1 MeV to (i.e. of the order of 0.5-0.3%). Gluonic corrections effectively 
change the pole mass m& ~ 4.7GeV to the running mass at scale mz '■ m^ — > rh\,{mz)- 
We calculate m& in terms of pole mass m& , a s (mz) and a s {m^) using standard two loops 



equation in MS scheme (see e.g. ref. [22]). 

For Z — > cc decay the running mass rh c (mz) is expected to be of the order of 0.5 GeV 
and corresponding contribution to T c is of the order of 0.05 MeV. Actually we take into 
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account this tiny effect in LEPTOP because it is incfuded in the other computer codes 
(ZFfTTER). 

In connection with T c , let us note that the term ^{t) given by eq. ( |1.20|) contains inter- 



ference terms which appear at the order (^) 2 , hence its negative sign. These terms have 
three types of final states: one quark pair, one quark pair plus gluon, two quark pairs. 
The last final state constitutes about 5 % of the J 2 and at the present level of experimen- 
tal accuracy is negligibly small. However, in principle, it calls for special care, especially 
when the two pairs consist of quarks of different flavor, e.g. bbcc. Such mixed final states 
should be a subject of special negotiations between theorists and experimentalists. 

Asymmetries 

a) Forward-backward asymmetry into // channel is given by the equation 



where for light fermions 



44 = ^A e A f (1.22) 



A f = < % Af l7 * ■ t 1 - 23 ) 

(9A f ) 2 + (g Vf ) 2 



For 6-quark we take into account the effect of nonzero mass: 



2gAb9vb 

[v 2 9Ab 

where v is the 6-quark velocity 



Ab ~W A b + W-v 2 W Vb ] v > (L24) 



\ 



1 - A <}f z) . (1.25) 



m z 



As was already mentioned, it is impossible to separate a given quark channel from another 
one unambiguously, starting from the order (— ) 2 , just due to the possibility of creating 
an additional pair of "foreign" quarks. So we prefer not to consider terms of the order of 
(2i) 2 in asymmetries at all. In this approximation the ratio gv / 9a is not renormalized by 
the final state interaction. 

b) Longitudinal polarization of r-lepton 

P T = -A T (1.26) 

c) The relative difference of total cross-section at Z peak for left- and right-handed elec- 
trons colliding with unpolarized positron beam 

A LR = °JiZ^ = Ae (1.27) 
Other observables are defined in terms of T q and T/ by equations: 

a) Hadronic width (up to very small corrections) is given by the sum of five quark channels: 



T h = T u + T d + T c + T s + T b 



(1.28) 
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b) Total width: 

v z = r h + r e + r M + r T + 3T U (1.29) 

c) Peak cross section of e + e~-annihilation into hadrons: 



d) Ratios 



i27rr e r\ 

° h = W Z T% (L30) 



R c — , Rb — ) Ri — ■ (1-31; 
J- h J- ft J- ; 



1.3 The a— Born approximation 

If we write the expressions for observables in the Born approximation in terms of s and c 
we automatically include purely electromagnetic interaction. Thus we get: 

(m w /m z ) B = c (1.32) 

(g Af ) B = T y (1.33) 
(gvf/gAf) B = l-MQf\s 2 (1.34) 

Substituting (gv,A) B into equations for the widths and asymmetries from the previous 
section we calculate the whole set of observables in Z decay in the Born approximation 
for electroweak interaction, but take into account trivial pure electromagnetic and gluonic 
interaction in the final state and the effect of the running of e.-m. coupling constant. 
Rather luckily for LEPTOP these equations reproduce the precision experimental data 
for Z decay and for mw (Jan. 95) with the accuracy la — 1.5a ||, 0, @- (Note that the 
final state interaction factors (eqs. (|1.15|) - ( |1.16|) ) include a lot of very small terms. One 
can neglect mass terms, as well as the difference between the vector and axial channels eq. 
( |1.16|) . Such simplified Born approximation will still reproduce the data with practically 



the same accuracy). 

It is important that the Born approximation for "low" -energy observables (parameters 
of Z boson decay and mw) is the same for very different models, hence the present 
experimental confirmation of the a-Born approximation is not sufficient to prove the 
validity of the MSM or to choose between different generalizations. 

1.4 Electroweak loop corrections 

1.4.1 Hadron-free observables 

For hadron-free observables we write the result of one-loop e.-w. calculations in the form 
suggested in ref. [||: 

3c 

m w /m z = c + 32?rg2(c2 _ g2) aKn(t, h) 
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9 - = - 1/2 -64rf?^ (Ml) 



R = gvi/gAi = 1 - 4s 2 + A _ U2 aV R (t, h) 



(1.35) 

3 

47r(c 2 — s 2 
g v = l/2+ 22 aV„(t,h) 

where t = m 2 /m 2 z , h = m 2 H /m 2 z and functions Vi(t, h) are normalized by the condition 

Vi(t,h)~t (1.36) 

at t > 1. 

Each function Vi is a sum of five terms 

K(t, fc) = t + Tt(t) + Hi(h) +Q + SVi(t, h) (1.37) 

The functions t + Tj(£) are due to (t, b) doublet contribution to self-energies of the vec- 
tor bosons, Hi{h) is due to W ± ,Z and H loops, the constants Cj include light fermion 
contribution both to self energies, vertex and box diagrams. 

To give explicit expressions for Tj(t) and H{(h) it is convenient to introduce three auxiliary 
functions F t (t), Fh{h) and F' h {h) (see subsection 4.3 for their expressions). 

The equations for Tj(i) and Hi(h) have the form ||: 

i = m 



m / x / 2 8 n v , 4 32 
rmW = ( 3 _ 9 " ' lni_ 3 + "9" a + 

1 2 _ 16 _ 2 _ 32 

3 c 4 3 c 2 L 3 9 3 9 J W ' V ' 

rr /, \ h , , c 2 h , h s 2 , 8 n 

# m (/i) = m/i + ^ln— r/i s 2 + 

v ; /i-l h-c 2 c 2 18c 2 3 

,/i 2 Ah 4,_ / ,, 
+( Y -y + 3)^)- 

+ (1.1205 - 2.595s 2 ) ; 
where 5s 2 = 0.23117 — s 2 (note the sign!). 

i = A 

TA(t) = I _ « > + 15 s , _ 1 - ««W 



3 9 27 4t - 1 
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■g s 4_ 8 2 1 4 2 l-2 tft(f ). 



TT / 1 \ c 2 , h 8h . . 

~TE h + [0-7751 + 1.075s 2 ] . 
18 



R 



T R (t) = hnt+^-^(l + llt)F t (t); 

Ah c 2 h 

H R (h) = — + -— hw + (1.40) 

RK 1 3 18 1 -c 2 /h c 2 K J 

+ (1 -th+ -h 2 )F h (h) + — h — In h+ 
v 3 9 9 ' K ' 1-h 

+ (1.3590 + 0.515s 2 ) 
i = v 
T v {t)=T A {t) 

H v (h) = H A (h) (1.41) 

The constants C; are rather complicated functions of sin 2 9 and we present their numerical 
values at s 2 = 0.23117 - 5s 2 

C m = -1.3500 + 4.135s 2 (1.42) 

C A = -2.2619 - 2.635s 2 (1.43) 

C R = -3.5041 - 5.725s 2 (1.44) 

C u = -1.1638 - 4.885s 2 (1.45) 



Functions 5V{(t, h) in eq. ( [1.371 ) are small corrections to They can be separated into 
five classes: 

1. Corrections due to iy-boson 5w& and t-quark 5 t a polarization of e.-m. vacuum 
are traditionally not included into the running of a(q 2 ). We also prefer to consider 
them together with electroweak corrections. This is especially reasonable because 
W contribution 5^a is gauge dependent. Here and for all other e.-w. corrections 
we use 't Hooft-Feynman gauge. The corrections 5w® and 5 t a were neglected in 
ref. || and were introduced in ref. fTj. 

16 1 

Stfmit, h) = Tis 4 -(5 w a + 5 t a) = -0.055 , (1.46) 

3 a 
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16 1 

SMt, h) = ns 2 c 2 -(8 w a + S t a) = -0.181 , (1.47) 

3 a 

6xV A (t,h) = 6iV v (t,h)=0, (1.48) 

where 

^ = — [(3 + 4c 2 ) (1 - V4C 2 - 1 arcsin —)--] = 0.0686 , (1.49) 
a 2n 2c 3 

5+a 4 r/ . „. . 1, 41 . 

"i" = "Si* 1 + 2( » F 'W - s> ^ "45^7 + - = - ' 00768 ' (L50) 



(Here and in the equations ( |1.51| ) - (1.69) we use m t = 175 GeV for numerical 
estimates.) 

2. Corrections of the order of aa s due to the gluon exchange in the quark e.-w. loops 
[ [Z4{ ] (see also |BJ). For the two generations of light quarks (q = u, d, s, c) it gives: 

SlV m (t,h) = 2 ■ [|(^^)(c 2 - S 2 )lnc 2 ] = (M5£))(_ .377) (1.51) 

S$V A {t,h) = %V v (t,h) = 2- [|(^^)(c 2 - S 2 + f /)] = (^^)(1.750) (1.52) 

%(t,/i) = (1.53) 

The result of calculation for third generation is rather complicated function of top- 
mass: 

%V m (t,h) = l(^p)){M x (l) + (1 - f + (± - ^ 2 )lnt 

-4(l-4)xtF 1 (^)-4^F 1 (0)} (1.54) 

C b C 

%V A (t,h) = %V v {t,h) = l(^l){M 1 (i) - 
+ (1 - ^ 2 ) 2 [^i(^) - ^i'(^)] + " ^ + ^ 4 ) - (0)} (i- 55 ) 

= - ^1(^-4^(0) + lint}, (1.56) 
3 7r 4r 3 4r 6 

where 

a s (m t ) = ® S [ M * I r - (1.57) 

sl " l + f^Mjln* 1 ; 

Note that 62 V%, unlike Vi themselves, do not depend on m#. The functions Vi(r), 
Ai(r) and i*i(x) have a very complicated form and were calculated in ref. |24J. For 
our purpose we can use their rather simple Taylor expansion for small values of their 
arguments (we have added cubic terms to the expansion presented in ref. [p4fQ : 
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V x {r) : 



: r 



r x 5, 2 328 , 1796 

■ 4C(3) - - + r 2 + r 3 

1 Sl ; 6 J 81 25 x 27 



+ . 



'1.56 



A 1 (r) = [-6C(3)-3C(2) + — ]+r[4C(3)-— ]+r 2 — + r 3 — + ... (1.59) 

H J i j sv;^ 4 j^ lsv; lg j -r 4Q5 ^ 7 x 25 x 27 v ; 



^19^ 1 95 

F 1 (x) = [--C(3) - -C(2) + -] + x[C(3) - -C(2) - -] + 

+ x 2 [-C(2) + -^—] + x 3 [—C(2) + —} + ... 
L 8 SV ; 3 x 64 J L 30 SV ; 72 J 



;i.60) 



where ((2) = tt 2 /6,C(3) = 1.2020569.... 

By summing contributions of ( |1.54|) -( |I.56| ) and using expansions (|1.58 )-( 1.60|) one 
obtains up to terms 0(l/t 3 ): 



%V m (t, h) = ( Q ^ m ^ )[-2.86t + 0.46 lnt- 1.540 



7T 



%V A (t,h) = 6 t 2 V v (t,h) 



,a s [m t 



71 



)[-2.86t + 0.493 



SlV R (t,h) 



,a s [m t 



7T 



)[-2.86t + 0.22 hit - 1.513 



0.68 


0.21 


_ a s (m t ) 


t 


t 2 J 


7T 


0.19 


0.05 


_ a s (m t ) 


t 


t 2 J 


7T 


0.42 


0.08 


_ a s (m t ) 


t 


t 2 J 


IT 



-11.67) 
(1.61) 

-10.10) 
(1.62) 

-11.88) 
(1.63) 

As these formulas are valid for m t > mz, in order to go to the region m t < mz we 
either put S^Vi = or use massless limit in which S^Vi = \8%Vi- m an y case this 
region gives tiny contribution into global fit. 

3. Corrections of the order of aa 2 were calculated for leading term aa 2 t only |25j 

5 3 Vi{t, h) ~ -(2.38 - 0A8N f )a 2 s (m t )t ~ -lA8a 2 s (m t )t = -0.07 (1.64) 
for Nf = 5 light flavors. (For numerical estimate we use a s (rriz) = 0.125.) 

4. The leading correction of the order of a 2 t 2 which originates from the second order 



Yukawa interaction was calculated in ref. ||T0|-P3| 

a . ,h 



5M(t, h) 



16ns 2 c 2 t 



AG)-* 



;i.65) 



Function A(j) is given in the table [T7T| JT^]-[|T^] for mn/m t < 4. We use expansion 
from [JTTyi — ||X1|| for mujm t > 4: 



A(- 



49 



27 



1 



+ 7r 2 H lnr + -ln 2 r + -r(2 - 12tt 2 + 121nr - 271n 2 r) + 

4 2 2 3 
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+ -^r 2 (1613 - 240vr 2 - 1500 In r - 720 In 2 r), (1.66) 

where r = t/h. For m t = 175 GeV and m# = 300 GeV one has A = 8.9 and 
S 4 Vi(t, h) = -0.11. 

5. In the second order in e.-w. interactions there appears quadratic dependence on the 
higgs mass |26| . 

6c ,m%r s 0.747 

^=2S ( 4 )X ~ =0 - 0011 (L67) 

n m 2 rr 1 1QQ 

5 5 Va = 5 S V V = ^(Hf) x -f^ = 0.0057 (1.68) 
247T m z s z 

S 5 V R = --^(^4)^f^ x 0.973 = -0.0032 . (1.69) 
247T m z s z c z 

(For numerical estimates we use m# = 300 GeV.) 
1.4.2 Z boson decay into hadrons 

For the partial width of the Z decay into a pair of quarks qq (q = u, d, s, c, b) we use the 
equation ( |1.14j ) 

T g = V2[g\ q R,A q + gy q Rvq]To , 
where e.-w. radiative corrections are included into gv q and gAq 

9M = T 3,[l + 32^ V ^(*» h )} ( L7 °) 
3IO I 

g Vq /g Aq = 1 - A\Q q \s 2 + aV^fo h) . (1.71) 

The functions VA q {t,h) and V Rq (t,h) in the one-loop e.-w- approximation are related to 
the functions V^t, /i) and Vfi(t, /i) from leptonic decays ||: 

1287TS 3 C 3 

V*,(t, /i) = Vk c (t, /i) = V A (t, h) + — — (F A i + F Au ) (1.72) 

da 

1287TS 3 C 3 

VA d (t, /i) = V As (t, h) = V A (t, h) + — — (Fai - F Ad ) (1.73) 

3a 

V Ru (t, h) = V Rc (t, h) = V R (t, h) + 16nsc [ c2 _ - fi2 ) x (1.74) 

3a 

x [F V i - (1 - 4s 2 )F Ai + ^(-(1 - ^s 2 )F Au + F Vu )\ 

\6tTSc(c 2 — S 2 ) 

V M {t, h) = V Rs (t, h) = V R (t, h) + \_ - } x (1.75) 

da 

x [F vl - (1 - 4s 2 )i^ + 3((1 - is 2 )F Ad - F Vd )} 
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where: 



F M = 


— (3.0088 + 16.45s 2 ) , 


(1.76) 


Fyi = 


-^(3.1868 + 14.95s 2 ) , 

4:71 


(1.77) 


Fau = 


--^(2.6792 + 14.75s 2 ) , 

47T 


(1.78) 


Fvu = 


- — (2.7319 + 14.25s 2 ) , 
An 


(1.79) 


FAd = 


— (2.2212 + 13.55s 2 ) , 

47T 


(1.80) 


Fyd = 


— (2.2278 + 13.55s 2 ) . 

All 


(1.81) 



The five digit accuracy of the above numbers is purely arithmetical. The physical uncer- 
tainties (caused by the neglect of the higher loop corrections) are at the level of the third 
digit. 

The difference V,i q — Vi (where % — A, R) is due to different e.-w. corrections to the vertices 
Zqq and Zll. 

The oblique corrections of the order of & s , 6> 2 s t and at 2 to the VAqiVRq) are the same as 
in the case of Va and Vr. But for Z-boson decay into pair qq there are also additional 
a s corrections to the vertices that have not been calculated yet. That brings additional 
uncertainty into the theoretical accuracy. 

For Z — > bb decays we have to take into account corrections to the Z — > bb vertex which 
depend on t ffifl, p§ 



3s 2 c 2 



V Ab (t, h) = V Ad (t, h) - — — — (0(f) + 6<f>(t)), (1.82) 
Vnit, h) = V Rd (t, h) - ^]f_~^ W) + W)) (1-83) 



where for <fi(t) we use the following expansion 

3 -2s 2 f o t 



0(t) = ^^{t + c 2 [2.88/n- -6.716- 



+ i(8.368c 2 ln 4r - 3.408c 2 ) + ^(9.126c 4 ln 4r + 2.26c 4 ) + (1.84) 
t c l t l c A 

+ 4r(4.043c 6 In 4t + 7.41c 6 ) + ... 
t 6 c l 

and for 5<f>(t) we use the leading approximation calculated in ref. [T2I| and JT0||- [|1~3]] . 

, , 3 — 2s 2 f 7r 2 .aAmf) . 1 .a, , e>\,hA 



1.4. Electroweak loop corrections 
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where the function rf® is given by the table |1.1| |TIJ-[]TB| for m#/m t < 4. For m#/m t > 4 
we use the expansion |I0|-O]: 



r b (2) (^) = ^[311 + 247r 2 + 2821nr + 901n 2 r-4r(40 + 67r 2 + 151nr + 181n 2 r) + 



+ r 2 (24209 - 6000tt 2 - 45420 lnr - 18000 In 2 r)L 

100 v n - 

where r = t/h. For m t = 175 GeV and m# = 300 GeV 

r? } = 1.245 . 



;i.86) 



Asymmetries are calculated by using eqs. Ql.22 ) - ( 1.27 ) with the loop corrected values 
of g A and g v . 



1.4.3 Auxiliary functions F t and 

The functions F t (t) and Fh(h) are used in eqs. (|1.38| ) - (|1.40|) . 



2(1 - ^/At^l 



arcsm 



Ft(t) 



At > 1 



2(1 - JT^At In , At < 1 



;i.87) 



F h {h) 



l + (^-|)ln h + hJl-lHJ%-l + J%) , h>A 



1 + 



h 

h-1 2' 



In h — hJ | — 1 arctan , h < A 



-1 + ^iln /, + ( 3 -/ i )J 7 ^ln( v /f-l + 



1 + ^f 1 In ft, + (3 - h)J^ arctan( 



ft < 4 



1.89) 
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Table 1.1: Functions A(^) and T ^{^) from @-P 



TTLf 


TTLf 


T (2)/S] 

fTlf. 


.00 


.739 


5.710 


.10 


1.821 


4.671 


.20 


2.704 


3.901 


.30 


3.462 


3.304 


.40 


4.127 


2.834 


.50 


4.720 


2.461 


.60 


5.254 


2.163 


.70 


5.737 


1.924 


.80 


6.179 


1.735 


.90 


6.583 


1.586 


1.00 


6.956 


1.470 


1.10 


7.299 


1.382 


1.20 


7.617 


1.317 


1.30 


7.912 


1.272 


1.40 


8.186 


1.245 


1.50 


8.441 


1.232 


1.60 


8.679 


1.232 


1.70 


8.902 


1.243 


1.80 


9.109 


1.264 


1.90 


9.303 


1.293 


2.00 


9.485 


1.330 


2.10 


9.655 


1.373 


2.20 


9.815 


1.421 


2.30 


9.964 


1.475 


2.40 


10.104 


1.533 


2.50 


10.235 


1.595 


2.60 


10.358 


1.661 


2.70 


10.473 


1.730 


2.80 


10.581 


1.801 


2.90 


10.683 


1.875 


3.00 


10.777 


1.951 


3.10 


10.866 


2.029 


3.20 


10.949 


2.109 


3.30 


11.026 


2.190 


3.40 


11.098 


2.272 


3.50 


11.165 


2.356 


3.60 


11.228 


2.441 


3.70 


11.286 


2.526 


3.80 


11.340 


2.613 


3.90 


11.390 


2.700 



Chapter 2: LEPTOP THEORETICAL 
UNCERTAINTIES 



2.1 General remarks 

With the increasing experimental accuracy of the precision measurements at LEP and 
SLC the question of the accuracy of the theoretical predictions becomes one of great 
importance. There are several major sources of uncertainties: 

1. The degree of accuracy with which any given observable is extracted from the exper- 
imental data. This involves theoretical models and computer simulations (Monte Carlo 
etc.) which deal with purely electromagnetic and strong interactions, for instance, in ex- 
tracting raw from the pp-collider data, or Z-observables (in particular for separate quark 
channels) from e + e~-data. We call them "extraction" uncertainties. 

2. Uncertainties in the input parameters (a, m z , G^, m b , m T etc.), so called "parametric" 
uncertaities. 

3. The degree of accuracy with which the theoretical expression for a given electroweak 
observable is derived within the MSM or some of its generalizations ("theoretical" uncer- 
tainties). 

In this paper we shall mainly deal with the third kind of uncertainties (in the framework 
of the MSM), although for some of observables, such as R b or A h FB or Qfb, they may not 
be fully separable from those of the first kind. One also has to keep in mind that some 
of the physical quantities we shall discuss such as Th,T e ,gA, are, strictly speaking, not 
primary, but derived observables. 

There are several types of observables with specific uncertainties for each type: 

1. " hadron-free" observables: raw { m w/ m z), T/j 9a, Qv/Qa- 

2. Observables of hadronic decays of Z boson. 

3. Observables of the low-energy weak processes: z/e-scattering, z/iV-scattering, atomic 
parity violation parameters. 

Let us finish this introductory section with the estimates of the "parametric" uncertainties 
for the Z-decay parameters and the W boson mass. Among the three basic input param- 
eters a , G^ and ra z the largest uncertainty comes from a = a(ra z ) : a -1 = 128.87(12). 
The relative uncertainty in ra\ is an order of magnitude smaller (mz = 91.1887(44)); that 
in G M is 50 times smaller (G^ = 1.16639(2) • 10~ 5 GeV~ 2 ). Induced by the uncertainty in 
a the absolute uncertainty in m w /ra z is close to 2 • 10~ 4 , that in g v / 1 9a is close to 1 • 1CT 3 , 
in T h it equals 0.8 MeV and in a h it is 0.001 nb. Finally, in T b it is 0.15 MeV. All these 
uncertainties come mainly from the uncertainty of the a-Born approximation due to that 
in the parameter s 2 = sin 2 ^ = 0.2312(3). 

Let us stress that relative parametric uncertainty in Qv/9a is quite large: about 1.5%. The 
experimental accuracy with which gv / 9a is measured is becoming close to the uncertainty 
due to 5a. Improvement of the accuracy of a calls for new measurement of e + e _ — > 
hadrons cross-section below J/ip resonanse. Measurement of (g — 2) M with high accuracy 
will also help. 
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2.2 Uncertainties in Vi 



Theoretical uncertainties come from as yet not calculated Feynman diagrams or not calcu- 
lated terms in a given diagram. The smallest calculated terms in the W and Z selfenergies 
are the virtual higgs corrections of the order a\yt 2 [jTH] - [ ]~3 1 , and the two gluon correction 



aw&it HH to the top quark loop. These corrections produce universal shifts 8Vi in the 
functions Vi which describe radiative corrections in the approach developed in 0: 



5V! 



a A(m h /m t )t 2 
7T 16s 2 c 2 



(2.1) 



-1.2a 2 (m t )t 



(2.2) 



Here and in what follows we denote by 5 the corrections, while by A - the uncertainties. 

According to f29j, still not completely calculated correction to ( |2.1|) , of the order of at, 
can be numerically close to ( |2.1| ). The same, according to the literature, is valid for un- 
certainties in eq. (|2.2|) . In order to have correct asymptotic behaviour of the uncertainties 
at t ^> 1 we assume that uncertainties correslponding to expressions ( |2.1| ) and (2]2) are 
derived by multiplying both of them by 2/t. For m# = 300 GeV , m t = 175 GeV 
and a s (m z ) = 0.125 our estimates are AV? = ±0.06 ,AVf s = ±0.03; corresponding 
uncertainties in the observables are presented in the Table 1. 



2.3 Uncertainties in hadronic Z decays 

In hadronic decays we separate decays into light quarks qq (where q = u,d, s, c) and the 
decay into bb. Concerning T q the next not yet calculated Feynman diagrams are the gluon 
corrections to the Zqq electroweak vertex "triangle" involving W or Z. The electroweak 
corrections to the vertices in the Z — > qq decays are given by the following equation: 

5T q = 24scT 2T«{(l-4\Q f \s 2 )F Vq + F Aq } . (2.3) 

Substituting numbers in eq. ( |2.3j ) we get: 

Ar„ = Ar c = -1.9 MeV , (2.4) 

AT d = AT S = -2.0 MeV . (2.5) 

Taking into account 4 light quark flavors and multiplying the results ( |2.4|) and (|2.5| ) by 
o. s {mz)/'^ we get the following estimate of the uncertainty in Z width into light quarks 
due to the gluonic perturbation of Zqq vertex: 

4 

J2 AT T S = ±0 - 3 MeV ( 2 - 6 ) 
i 



For the Z — >• bb decay the leading term ~ ay/& s t [p8 |, and the potentially next to the lead- 
ing term ~ a^a s lnt [30] come from virtual gluons in the ttW vertex triangle. Calculated 



2.4. Procedure to estimate total theoretical accuracy of LEPTOP 
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28fl the leading a s correction to electroweak vertex lead to the following correction to 



m 

rv 

-,2\ ^2 



Arp* = -r r o n,2 2 (-— )-^* = °- 7 MeV > ( 2 - 7 ) 



a (3 - 2s 2 ) ir 2 a s (m t 

~ — 2 2 1 5" J 

where for numerical estimate we substitute m t = 175 GeV. 

Calculated in [M d s correction to In t term appears to be much smaller: 



5T m s mt = _<* r (3 l s y i^?A lnt „ Q.002 MeV (2.8) 

7T 16S J C 2 81 7T 

If we suppose that uncalculated terms are of the order of the last calculated one we can 
safely neglect this type of uncertainties. 

Uncertainty in also shows itself in Ft, value. However, it is comparatively small. Even 
if we assume that is known with accuracy ±300 MeV, the shift in r& is ±0.17 MeV. 
There is a specific to Z — > bb decay virtual higgs gaugeless correction to the vertex 5(f>(t, h) 
consisting of terms of the order a^t 2 , a^t In t , a^t and so on. Only the first of them 
have been calculated [I(|-[|l3| : 

If we consider this leading term as an upper bound of the "higgs theoretical uncertainty" 
in A(p than we get: 



, , a . a 3 — 2s 2 at 2 ,■■>, 

ATI = --IVtf = --r ——rl» = 0.02 MeV , (2.10) 



where we substitute m t = 175 GeV , trh = 300 GeV. This uncertainty in r& is much 
smaller than those considered above and could be neglected. However, correction of the 
order of a 2 (m t )t to (j)(t), 5(j) as , may lead to larger correction to IV 

To estimate its value let us suppose that it's ratio to the calculated correction ~ 6t s t is 
equal to the ratio of the calculated corrections ~ a 2 t and ~ a s t to functions Vf 



" & £ ^ i )(_2. 86 )i 
3 - 2s 2 1.37?r 



2s 2 c 2 



a 2 s (m t )t. (2.12) 



Taking this estimate as a source of uncertainty in T b , for m t = 175 GeV we get: AI^ S 
0.1 MeV, which dominates over fl2.8|) and (p,10|). 



2.4 Procedure to estimate total theoretical accuracy of LEP- 
TOP 

In order to estimate the total theoretical uncertainties for the observables we implement 
in LEPTOP the procedure proposed by D.Bardin and G.Passarino |HJ]. We choose several 
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options to the preferred formulas used in LEPTOP and make variations of these formulas. 
Usually each option consists in the addition of an extra term corresponding to a rough 
guess of the value of the uncalculated higher order terms. We then make all possible 
combinations of these options, i.e. 2 n in total, where n is the number of options. Among 
all these 2 n combinations we locate those yielding the minimum and the maximum values 
of the observables and took as the estimate of the theoretical errors their deviations from 
the central values. 

We consider the following options: 

Options 1,2 Correction SVf given by eq. Q2.1| ) multiplied by 2/t is added to (option 
1) (or subtracted from (option 2)) the three functions Va , V m and Vr. 

2 

Options 3,4 The same as in options 1,2 with SV"" from eq. (|2.2| ). 
Options 5,6 The expression for T q is modified by adding (option 5) or subtracting 
(option 6) YliAT q = 0.3 MeV in order to take into account gluon cor- 



Option 7 



rections to electroweak Zqq triangle vertices. 
Correction 50°^ eq. (|2.12|) is added to function 4>{t). 



2.4. Procedure to estimate total theoretical accuracy of LEPTOP 
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Table 2.1: The effect of the working options of LEPTOP on theoretical errors. The first 
two lines indicate parametric uncertainties caused by 5s 2 = 0.0003 (which is equivalent to 
5a" 1 = 0.12) and by Suit, = 0.3GeV. The next six lines refer to the intrinsic theoretical 
uncertainties. This table was calculated with m t =175 GeV, Mh=300 GeV and d s =0.125. 





lib yy 


r. 


Sill 2 l r r 

bill U e ff 




IV 

L Z 


1 h 




L 


± " 


R x 




MeV 


MeV 






MeV 


MeV 




MeV 


10 5 


10 5 


Ss 2 


16 


.015 


.00031 


.0014 


.8 


.8 


.0055 


.15 


1.1 


1.8 










.0023 


.2 


.2 


.0029 


.24 


11.0 


2.4 


AVf 


9 


.015 


.00005 


.0002 


.5 


.4 


.0009 


.08 


.2 


.3 




5 


.008 


.00003 


.0001 


.3 


.2 


.0005 


.04 


.1 


.2 


AT q 








.0029 


.3 


.3 


.0037 




3.8 


1.5 










.0010 


.1 


.1 


.0012 


.10 


4.6 


1.0 


total 


+13 


+.023 


+.00008 


+.0032 


+1.2 


+ 1.0 


+.0063 


+.23 


+8.7 


+ 1.9 




-13 


-.023 


-.00008 


-.0042 


-1.1 


-.9 


-.0051 


-.13 


-4.1 


-2.9 
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3.1 Introduction 

LEPTOP code allows the calculations of LEP-I electroweak observables in the framework 
of the Minimal Standard Model and fitting the experimental data for determination of 
m t , M H and a s (Mz)- LEPTOP package has been written in the Fortran language and 
installed on IBM VM and HP-SUN UNIX operating systems. It uses the subroutines 
from the CERN program library f32j. The code is maintained with the CMZ [[33| source 
code management system. The Fortran code is available either in the form of a compilable 
source file leptop.f or in the CMZ ASCII readable file leptop.car. In the following sections 
we describe simple useage of the LEPTOP package. The user interface part of LEPTOP 
was written in the object oriented style, which allows a simple use of it without deep 
knowledge of the Fortran language. The data encapsulation principle and hiding the 
internal variables are applied. This significantly reduces the amount of information the 
user should care about improves the robustness of the program. 

In order to obtain the maximum accuracy we used whenever possible double-precision 
variables in the internal calculations. However the user interface was written in a single 
precision in order to simplify the use of the program, make easy interface with other 
standard packages and avoid user's bugs typical for the mixture of single and double 
precision variables. The 7 digit precision (for 32 bit machines) in the input and output 
is adequate for the present experimental and theoretical errors. The units used in the 
program are GeV and nanobarns, unless explicitly stated. 

This document is a Reference manual and User's Guide. Whener possible we tried 
to follow cernman style and notations . An example of a simple program is given in 
the section 3.8. 

This report, including the writeup and Fortran code, is available under WWW [3lJ on 
http : / /cppm . in2p3 . f r/leptop/ intro_leptop . html. 



3.2 Initialisation 



CALL LTINIT (ISTAT) 

Action: The setting of the default values of all variables is done. LTINIT (0) should be 
called before any other LEPTOP routine. 

Input parameters: 

ISTAT Level of the initialisation. 

ISTAT=0 means complete initialisation from scratch. 

ISTAT=1 means recalculating constants after the eventual update of the data. 

This option is needed only for advanced users, because in the user 
interface it is done automatically. 



20 



3.3. Input and output 



21 



3.3 Input and output 

The output of the program is directed into the unit 6. The unit 19 is reserved for the 
internal scratch output. The amount of output is controlled via PRNT flags (see LTFLAG 
routine) . 

3.4 Setting of flags and options 

CALL LTFLAG (KEY, ID) 

Action: Setting on and off the flags according to the character variable KEY, identifying 
the group of flags, and integer index ID, identifying specific flag of that group. 
Positive values of ID sets the corresponding flag on, negative value of the ID sets it off. 
Zero value of the ID sets off all the flags of the corresponding group. By default all the 
flags and options are switched off. 

Input parameters: 

KEY Character variable identifying the group of flags. 

' PRNT ' Group of flags controlling the output level. 

' OPT ' Group of options modifying the prefered formulas of LEPTOP. See 
section 2.4 for details. 

ID Numerical identifier of the flag within the group. 
The meaning of different flags is given below: 
'PRNT' Allowed values for ID are: 

1 Low level printout for debug purpose in the process of calcula- 
tion. 

2 Detailed printout during the fit. It may produce very big out- 
puts. 

3 Detailed printout after the fit. Printout of the table with MSM 
and Born predictions for some observables in the style of the 
table in the paper |J. 

4 Printout during the initialisation of several default constants and 
experimental inputs. 

8 Printout during the setting of the variables and experimental 
data is done by the user interface routines LTPUT, LTFPUT, 
LTFCOR, LTFUSE etc. 

9 Printout of the calculated observables, setting constants and ex- 
perimental data is done by the user interface routines LTGET, 
LTFGET , LTFIT1 , LTFIT2 etc. 

'OPT' Allowed values for ID are: 

1 Option 1: functions Va, Vm, Vr are increased by the terms (2/t)8Vf . 

2 Option 2: same variation as option 1, but with negative sign. 

2 

3 Option 3: functions Va, Vm, Vr are increased by the terms (2/t)8V" a 

4 Option 4: same variation as option 3, but with negative sign. 
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5 Option 5: the quark widths of Z are increased in order to take 
into account estimated QCD corrections to Zqq triangle vertices. 

6 Option 6: same variation as option 5, but with negative sign. 

7 Option 7: function <f>(t) is increased by the term <50 a = 

8 Option 8: the b-quark width of Z is increased by 7 MeV. 
'MNUNIT' Changing the Fortran output unit of MINUIT to unit ID. The default 

is the the output on the scratch file on unit 19. Allowed values for ID 
are from 1 to 99. To restore the output on the scratch file one should 
assign ID to -19. 

' MNSAVE ' Changing the SAVE output unit of MINUIT to unit ID. The default 
unit is 7. 

' MNREAD ' Changing the Fortran input unit of MINUIT to unit ID. The default 
value of ID is -5, which means that MINUIT is working in Fortran- 
callable mode. Attention: using this flag switch MINUIT to the 
data-driven mode, when user should provide on unit ID the complete 
set of MINUIT commands. If needed, the user can provide his own 
version of the function FCNLB. 

' MNPRNT ' sets the print level of MINUIT. The default value of ID is -1, which 
means no output of MINUIT except from MINUIT SHOW command. 
Other possibles values of ID are in the range from to 3 according to 
MINUIT SET PRIntout command. 

' MNEPS ' informs MINUIT that the relative floating point arithmetic precision 
is 10 ID . The default precision is 1CT 10 . 

Remark: 

LTFLAG is the only routine that can be called before LTINIT(O). This allows to suppress 
the printout during the initialisation stage. 

Examples: 



CALL LTFLAG('OPT',7) 



This call sets on the option 7. 
CALL LTFLAG('PRNT',-4) 

This call sets off the printout of the fit results, if it has been set on previously. 



3.5 Setting the constants and experimental data 



3.5. Setting the constants and experimental data 
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CALL LTPUT (NAME, DATA) 

Action: Setting the constant denned by keyword NAME to the value defined by the variable 
DATA. 

Input parameters: 

NAME Character variable containing the name of the constant to be modified. 
'MT' Mass of the top quark m t . 

'MH' Mass of the Higgs boson M H . 

'MZ' Mass of the Z boson M z . 

'ALBAR' Electromagnetic coupling constant d(M z ). 
'ALSHAT' Strong coupling constant a s (M z ). 

' GFERMI ' Weak Fermi coupling constant determined from the muon lifetime 
measurements. 

DATA New value of the constant. 
Remark: 

If PRNT flag number 8 was set before the call to LTPUT the new value will be printed. 
Example: 

CALL LTPUT('MH',300.) 
This call will set the Higgs mass to 300 GeV. 



CALL LTFPUT (TYPE, NAME, DATA) 

Action: Setting the values of experimental data, including observables and their errors, 
for the use in the fit of experimental data with LEPTOP formulas. 

Input parameters: 

TYPE Character variable selecting the type of data to be set. 

'VALUE' The central value of the experimental data is set. 
'ERROR' The error on the experimental data is set. 

'MARSEILLE' Experimental data presented at Marseille-93 conference are loaded 

in case NAME= ' ALL ' . 
'M0RI0ND94'Same, but for Moriond-94 conference. 
'GLASGOW Same, but for Glasgow-94 conference. 
'M0RI0ND95' Same, but for Moriond-95 conference. 
NAME Character variable identifying the experimental observable. The following ob- 
servables may be used in the fit: 
' GZ ' Total width of Z boson T z . 

'SIGH' The hadronic pole cross section a° h = {\2n / M 2 z )(T e T h /T 2 z ). 
'RL' Ratio Ri = T h /Ti of hadronic to leptonic partial width of Z. 

'AFBL' Forward-backward lepton asymmetry at Z pole. 
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' S2NUN ' 



'RC 
'ALL' 



'MTV 



'ALR' 



'QFB' 



'RB' 



'AFBC 



' ATAU ' 
' AETAU ' 
' AFBB ' 



'MWMZ' 



Ratio Rb = Tb/Th- 

Ratio M w /M z of W-boson to Z-boson masses. 

A T from r polarization. 

A e from r polarization asymmetry. 

Forward-backward b-quark asymmetry at Z pole. 

Forward-backward c-quark asymmetry at Z pole. 

Measurement of forward-backward quark asymmetry translated in a 

value for sin 2 Q e ff- 

Measurement of the ratio of neutral to charged current cross-sections 
translated into sin 2 6w- 

Left-right polarization asymmetry A LR , measured at SLAC by SLD. 

Top quark mass m t measured by CDF and DO at Fermilab. 

Ratio R c = T c /Th of charm to hadron partial widths. 

used to set all observables from the data set defined by the keyword 

TYPE. 



DATA Numerical value of the experimental data. 
Remark: 

By default Glasgow-94 experimental data are loaded during the initialisation with LTINIT(O) . 
Examples: 

CALL LTFPUT('VALUE','GZ',2.4974) 

This call sets the experimental value of Y z to 2.4974 GeV. 
CALL LTFPUT('ERROR','SIGH',0.12) 

This call sets the experimental value of the error on the peak hadron cross section to 0.12 



CALL LTFPUT('GLASGOW', 'ALL', dummy) 
This call sets the input of all experimental data as presented at Glasgow-94 conference. 

CALL LTFCOR (NAME1 , NAME2 , DATA) 

Action: Set the values of correlation coefficients on the experimental errors, for the use 
in the fit of experimental data with LEPTOP formulas. 

Input parameters: 

NAME1 Character variable identifying the first experimental observable, as defined in the 
LTFPUT description. 

NAME2 Character variable identifying the second experimental observable, as defined in 

the LTFPUT description. 
DATA Numerical value of the correlation coefficient between the two experimental data. 



nb. 



3.6. Access to the constants, observables and experimental data 



25 



Example: 

CALL LTFCOR('GZ','SIGH',-0.110) 
This call set the correlation coefficient between the experimental errors of T z and cr° equal 

t0 Pr z< = -°- U - 



3.6 Access to the constants, observables and experimental data 

CALL LTGET (NAME, DATA*) 

Action: Extracting the constant or the theoretical prediction for the observable defined 
by keyword NAME to the variable DATA. 

Input parameters: 

NAME Character variable containing the name of the constant or theoretical value of 
the observable to be extracted. 

'MT' Mass of the top quark m t . 

'MH' Mass of the Higgs boson M H . 

'MZ' Mass of the Z boson M z . 

'ALBAR' Electromagnetic coupling constant d(M z ). 

'ALSHAT' Strong coupling constant a s (M z ). 

' GFERMI ' Weak Fermi coupling constant G? M defined by muon lifetime measure- 
ments. 

' GV ' Vector coupling of the electron. 

' GA ' Vector axial coupling of the electron. 

'SIN2E' The effective electro-weak mixing angle sin 2 # e // for leptons. 

'SIN2B' The effective electro-weak mixing angle sin 2 e // for b quarks. 

' MW ' Mass of the W boson M w . 

'GNU' Partial width of Z boson into vv. 

'GE' Partial width T e of Z into e + e~. 

' GMUON ' Partial width of Z into 

' GTAU ' Partial width T T of Z into t + t~. 

'GU' Partial width r„ of Z into uu quarks. 

' GD ' Partial width T d of Z into dd quarks. 

' GS ' Partial width T s of Z into ss quarks. 

' GC ' Partial width T c of Z into cc quarks. 

' GB ' Partial width r& of Z into bb quarks. 

'GINV Invisible width of Z boson T inv . 

'GH' Hadron width of Z boson T h . 

' GZ ' Total width of Z boson T z . 

'RL' Ratio Ri = T h /Ti of hadronic to leptonic partial width of Z. 
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'RC' 


Ratio R c = T c /Th of charm to hadron partial widths. 


'RB' 


Ratio R b = T b /T h . 


' AFBL ' 


Forward backward lepton asymmetry A l FB . 


' AFBC ' 


Forward-backward c-quark asymmetry at Z pole. 


' AFBB ' 


Forward-backward b-quark asymmetry at Z pole. 


'ALR' 


Left right asymmetry A LR . 


'SIGH' 


The hadronic pole cross section cr° = (127r/M^)(r e r^/r^) 


Output parameter: 



DATA numerical value of the constant or the result of the calculation for the prediction. 
Remark: 

- The extracted value may depend on the change of the constants by the previous 
call to LTPUT. So to get a consistent set of observables one should make all calls to 
LTGET after the last call of LTPUT. 

- If PRNT flag number 9 was set before the call to LTGET the value of the observable 
will be printed. 

Example: 

CALL LTGET ('MW'.AMW) 
After this call the variable AMW will contain the W boson mass in GeV. 



Action: Extracting the values of experimental data, including observables and their 
errors. It also can extract the values of the fitted parameters and their errors after the 
fit. 

Input parameters: 

TYPE Character variable selecting the type of data to be extracted. 

'VALUE' The central value of the experimental data is given. 

'ERROR' The error on the experimental data is given. 

'FIT_RESULT' The central value of the fitted parameter after the last call to 
LTFIT1 or LTFIT2. 

'FIT_ERR0R' The corresponding parabolic error of the fitted parameter. 
'FIT_ERR+' The corresponding positive error of the fitted parameter. 
'FIT_ERR-' The corresponding negative error. 

'FIT_GLB' The corresponding global correlation coefficient with other fitted 
parameter. 

NAME Character variable identifying the experimental observable as defined in the de- 
scription of LTFPUT for TYPE= ' VALUE ' or TYPE= ' ERROR ' . If NAME= ' ALL ' and PRNT 
flag 9 is set on, the table of all experimental data is printed. If TYPE='FIT_*' 
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the allowed values of NAME are: MT, MH and ALSHAT, corresponding to the fitted 
parameters. 
Output parameters: 

DATA Numerical value of the experimental data or fitted parameter. 
Remark: 

- If PRNT flag number 9 was set before the call to LTFGET the value of the observable 
will be printed. 

Examples: 

CALL LTFGET ( 'VALUE' , 'MWMZ ' ,RMWMZ) 

After this call the variable RMWMZ will contain the experimental ratio of W boson to Z 
boson masses M w /M z . 

CALL LTFGET('ERROR','SIGH',ESIGH) 

After this call the variable ESIGH will contain the experimental value of the error on the 
peak hadron cross section a° in nb. 

CALL LTFLAG('PRNT',9) 

CALL LTFGET (' ',' ALL', DUMMY) 

After these calls the table with all experimental data is printed. 

CALL LTFITl('MH',amh,emh,chi2) 

CALL LTFGET('FIT'ERR+','MH',emhpos) 

CALL LTFGET('FIT'ERR-','MH',emhneg) 

After these calls the variables emhpos and emhneg will contain the asymmetric errors on 
Higgs mass, which might be significantly different from the parabolic error emh due to the 
non-gaussian distributions. 



3.7 Fitting experimental data with LEPTOP formulas 



The fitting in LEPTOP is done using MINUIT package ]36|]. The fitted parameter values 
correspond to the minimum of the \ 2 function: 



x 2 = E 



^ {Ei - Ti (m t , M H ,a s )) {E j - Tj (m t , M H , a s )) 



S, 



1 



-1 



Vij = pijSEidEj 
where the following notations were used: 



(3.1) 

(3.2) 
(3.3) 
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n Number of experimental observables used in the fit. 

Ei Experimental observable, which usually comes from the compilation of results 

of all LEP detectors. 

Ti Theoretical prediction of LEPTOP as a function of parameters m t) M H and a s . 

Vij Covariance matrix of errors. 

5Ei The experimental error on the observable i. 

Pij The correlation coefficient between the experimental errors of the observable i 

and observable j . 



The initial values of parameters are used from the settings done by the routine LTPUT. 
So the user can check the sensitivity of the fit to the initial values of the parameters. It 
is the user's responsibility to provide reasonable values of the initial parameters, as the 
minimisation procedure may diverge and move the parameters to the unphysical region, 
where the program might be not protected. 



Action: Instructing LEPTOP to use or not to use the particular experimental data in 
the fit. By default all experimental data listed in LTFPUT description are used, except 
'MT\ 

Input parameters: 

COMMAND Character variable selecting the command. 

'USE' Command to use this data in the fit. 

'NOUSE' Command to exclude this data from the fit. 

NAME Character variable identifying the experimental observable as defined in the 

description of LTFPUT. If NAME= ' ALL ' all experimental data are affected. 

Examples: 

CALL LTFUSE('USEVMT') 

Add the experimental value of m t from CDF to the fit. 
CALL LTFUSE('NOUSEVQFB') 

Do not use the experimental value of Qfb in the fit. 

CALL LTFIT1 (PARAM, PARI* , EPAR1* , CHI2*) 
Action: 

Performing one parameter fit of the experimental data. The actual parameter is selected 
by the character constant PARAM from m t , M H and a s . 

Input parameters: 

PARAM Character variable selecting the parameters. 
'MT' Select the parameter m t . 
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'MH ; Select the parameter M H . 

'ALS' Select the parameter a s . 
Output parameters: 
PARI Value of the fitted parameter. 

EPAR1 Estimate of the error of the fitted parameter. 
CHI2 Value of the x 2 of the fit. 

Example: 

CALL LTFIT1('MH',AMH,EMH,CHI2) 

After this call the variable AMH contains the fitted Higgs mass M H) the variable EMH 
contains the parabolic estimate of the error of the Higgs boson mass and CHI2 the x 2 of 
the fit. 



CALL LTFIT2 (PARAM , PARI* , PAR2* , EPAR1* , EPAR2* , R12* , CHI2*) 
Action: 

Performing two parameter fit of the experimental data. The actual parameters are selected 
by the character constant PARAM from m t , M H and a s . 

Input parameters: 

PARAM Character variable selecting the parameters. 

'MT.ALS' Select the parameters m t and a s . 

' MT , MH ' Select the parameters m t and M H . 

'ALS^tF Select the parameters M H and a s . 
Output parameters: 
PARI Value of the first parameter. 

PAR2 Value of the second parameter. 

EPAR1 Value of the error of the first parameter. 
EPAR2 Value of the error of the second parameter. 
R12 Value of the correlation coefficient between two parameters. 

CHI2 Value of the x 2 of the fit. 

Example: 

CALL LTFIT2('MT,ALS',AMT,ALS,EMT,EALS,R12,CHI2) 

After this call the variable AMT contains the fitted top quark mass m t) the variable ALS 
contains the fitted strong coupling constant a s (M z ), the variable EMT contains the er- 
ror of the top quark mass, the variable EALS contains the error of the strong coupling 
constant, the variable R12 contains the correlation coefficient between the errors of two 
parameters and CHI2 contains the x 2 °f the fit. 
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3.8 Sample program 

program example 
call ltdemo 
end 

SUBROUTINE LTDEMO 

* template program to demonstrate LEPTOP package * 

* set all LEPTOP print flags off 

CALL LTFLAG ( ' PRNT ' ,0) 

* initialisation 

CALL LTINIT(O) 

* get some constants 

CALL LTGET ( ' MZ ' ,AMZ) 

CALL LTGET ( ' GFERMI ' , GFERMI) 

CALL LTGET ( ' ALBAR ' , ALBAR) 

* get more constants 

CALL LTGET ( ' MELE ' , ame 1 e ) 
CALL LTGET ( ' MMUO ' , ammuo ) 
CALL LTGET ( ' MTAU ' , amtau) 
CALL LTGET (' MS ' ,amstr) 
CALL LTGET ( ' MC ' , amchrm) 
CALL LTGET ( ' MB ' , ambot ) 

* modify some constants 

CALL LTPUTC'MT' ,175.) 
CALL LTPUTC'MH' ,300.) 
CALL LTPUT ( ' ALSHAT ',0.125) 

* get some physical quantities 

CALL LTGET ('GV ,gv) 
CALL LTGET ('GA' ,ga) 
CALL LTGET (' MW ,amw) 
CALL LTGET ('GNU' ,gnu) 
CALL LTGET (' GE ' ,ge) 
CALL LTGET ( ' GMUON ' , gmuon) 
CALL LTGET ( ' GTAU ' , gt au) 



* fit of mtop and alsbar 

CALL LTFIT2 ( ' MT , ALS ' , amt 1 , als 1 , eamt 1 , eals 1 , rho 1 , chi21 ) 
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* set data from Glasgow conference 

CALL LTFPUT ( ' GLASGOW ' , ' ALL ' , . ) 

* modify value of Z total width 

CALL LTFPUT ( ' VALUE ' , ' GZ ' , 2 . 4974) 

* fit of mtop and mh 

CALL LTFIT2 ( ' MT , MH ' , amt2 , amhl , eamt2 , eamhl , rho2 , chi22) 

* fit mt 

CALL LTFIT1 ( ' MT ' , amt3 , eamt3 , chi23) 

* fit mh 

CALL LTFIT1 ( ' MH ' , amh4 , eamh4 , chi24) 

* fit als 

CALL LTFIT1 ( ' ALS ' , als , eals , chi2) 



END 



! LTDEMO 
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APPENDIX A 

FLOWCHART OF LEPTOP 



Select the three most accurate observables: 

G^M z ,a{M z ) = a 



Define angle ( s = sin (9, c = cos 0) 
in terms 01 G^, M z and a: 

G, = (7r/V2)a/s 2 c 2 M 2 z 



Define the Born approximation for other electroweak 
observables in terms of G^, M z and 0. 



Introduce bare couplings (do, az , ctw ), masses (Mz , Mw , 
m qo (including m to )) and VEV rj in the framework of MSM. 



Express o?o, a%, Mz Q in terms of G M , M z and a in the one-loop 
approximation, using dimensional regularization (l/e,/i) 



Express one-loop corrections to all other electroweak observables 
in terms of a*o, az , Mz , mt, M H and hence in terms of G M , 

M z , a, m t , M H . Check cancellation of 1/e and \i. 



Introduce gluonic corrections into quark loops and QED and QCD 
final state interactions for hadronic decays 

(in terms of a,d s (mz)*,rnb(M z ),rnt). 



Compare the Born results and Born + one loop results with 
experimental data on Z - decays and M w . 



Make a global fit for three parameters m t , M H ,d s (M z ). 



Predict the central values of all electroweak observables 
and their uncertainties. 

* — d s (M z ) is the QCD coupling in MS - scheme. 
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List of LEPTOP papers § 

1. Electroweak radiative corrections and top quark mass, CERN-TH. 6053/91, March 
1991, TPI-MINN-91/14-I, ITEP-15/91. 

2. Parametrization of electroweak radiative corrections, ITEP-104/92, November 1992; 
ZhETF 103 (1993) 1489 (in Russian), May 1993, JETP 76 (1993) 725. 

3. On the electroweak one-loop corrections, CERN-TH.6538/92; ITEP-67/92, June 1992; 
Erratum August 1992; Nucl.Phys. B397 (1993) 35. 

4. On the interpretation of the CHARM II data, CERN-TH.6695/92, October 1992; 
Phys.Lett. B298 (1993) 453. 

5. Virtual gluons in the electroweak loops (with N.Nekrasov), CERN-TH. 6696/92, Octo- 
ber 1992; Yad.Fys. 57 (1994) 883. 

6. Do-it-yourself analysis of precision electroweak data, CERN-TH. 6715/92, November 
1992; Erratum March 1993; Phys.Lett. B299 (1993) 329; Erratum 304 (1993) 386. 

7. The isolines of electroweak radiative corrections and the confidence levels for the masses 
of the top and higgs (with V.Yurov), CERN-TH.6849/93, March 1993; Phys.Lett. B308 
(1993) 123. 

8. On the electroweak and gluonic corrections to the hadronic width of the Z boson, 
CERN-TH.6855/93, April 1993; Phys.Lett. B320 (1994) 388. 

9. Do present data provide evidence for electroweak corrections? CERN-TH. 6943/93, 
July 1993; ModernPhys.Lett. A8 (1993) 2529; Erratum 8 (1993) 3301. 

10. On the effective electric charge in the electroweak theory, CERN-TH. 7071/93, Novem- 
ber 1993; Phys.Lett. B324 (1994) 89. 

11. The values of mt and a s derived from non-observation of electroweak radiative cor- 
rections at LEP: global fit (with A.Rozanov and V.Yurov), CERN-TH. 7137/94, January 
1994; Phys.Lett. B331 (1994) 433. 

12. The Q 2 dependence of W and Z coupling constants in the interval < \q 2 \ < m 2 ?, 
CERN-TH.7153/94, January 1994; ModernPhys. Letters A9 (1994) 1489. 

13. First evidence for electroweak radiative corrections from the new precision data (with 
A.Rozanov), CERN-TH.7217/94, April 1994; ModernPhys. Letters A9 (1994) 2641. 

14. Do the present electroweak precision measurements leave room for extra generations? 
(with A.Rozanov and V.Yurov), CERN-TH. 7252/94, May 1994. 



V.Novikov, L.Okun, and M.Vysotsky are coauthors of all fourteen papers 
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